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In ponderomotive spectroscopy an amplitude-modulated optical standing wave is employed to 
probe Rydberg-atom transitions, utilizing a ponderomotive rather than a dipole-field interaction. 
Here, we engage nonlinearities in the modulation to drive dipole-forbidden transitions up to the 
fifth order. We reach transition frequencies approaching the sub-THz regime. We also demonstrate 
magic-wavelength conditions, which result in symmetric spectral lines with a Fourier-limited feature 
at the line center. Applicability to precision measurement is discussed. 


Measurements of atomic transition frequencies are the 
cornerstone of precision metrology, used in applications 
ranging from atomic clocks ^ to measuring gravitational 
redshifts [2] and the radius of the proton [3j. An im¬ 
portant metric of precision in these applications is the 
fractional frequency resolution, A\vjv, in which v is the 
measured frequency and A;/ is its uncertainty. In or¬ 
der to obtain the best it is desirable to increase 

V while decreasing A;/. In recently-developed pondero¬ 
motive spectroscopy [4], Rydberg atoms are trapped in 
a standing-wave laser field (optical lattice). Electronic 
transitions are driven by modulating the lattice-light in¬ 
tensity at the transition frequencies of interest. In this 
Letter, we employ nonlinearities intrinsic to this excita¬ 
tion process to increase v by driving atomic transitions at 
higher harmonics of the modulation frequency. Unlike in 
standard nonlinear spectroscopy, in ponderomotive spec¬ 
troscopy the intensity of the probe field does not need 
to be increased to drive these higher-order transitions, 
keeping light shifts low. We also identify magic transi¬ 
tions that minimize residual trap-induced shifts, in addi¬ 
tion to reducing the uncertainty Azz of the line center by 
about one order of magnitude. 

The ponderomotive interaction of a Rydberg electron 
with an optical field allows one to drive transitions when 
there is substantial spatial variation of the field inten¬ 
sity within the volume of the atom and when the lattice 
potential is modulated in time at a resonant transition 
frequency between states. The interaction is described 
by the minimal-coupling Hamiltonian (a.u.) [S] 

iL = A(f) • p-b A(f) • A(f)/2, (1) 

where p is the electron’s momentum operator and A the 
vector potential of the laser electric field. In this work, 
the A-A (ponderomotive) interaction drives atomic tran¬ 
sitions, rather than the usual A • p interaction, leading 
to flexible selection rules [6] and expanding the range of 
transitions available for precision measurement. 

In Fig. we show the experimental set-up. See [1] for 
details. A continuous-wave (c.w.) 1064-nm laser beam 
is split by a Mach-Zehnder interferometer into a low- 
power and higher-power beam. The low-power beam is 
amplitude-modulated via an EOM driven by a microwave 
signal with voltage amplitude and frequency Hi. This 



FIG. 1. Experimental set-up. A Mach-Zehnder interferom¬ 
eter splits and re-combines two 1064-nm c.w. beams. One 
beam is amplitude-modulated by a fiber-based electro-optic 
modulator (EOM), driven at frequency Hi (inset). The op¬ 
tical lattice is formed by retro-reflecting and focusing the re¬ 
combined beam into a ®®Rb magneto-optical trap (MOT). 
The lattice potential is modulated with period T. 


beam is coherently re-combined with the unmodulated 
(higher-power) beam to parametrically amplify the mod¬ 
ulation sidebands [the radical term in Eq. ([^]. We form 
a standing-wave optical lattice in the atom-field interac¬ 
tion region by retro-reflecting the lattice beam. Using 
a lattice inversion technique 0, we laser-excite Rydberg 
atoms such that their center-of-mass locations are either 
at lattice intensity minima or maxima. Intensity modu¬ 
lation of the lattice then results in a time-periodic atom- 
field ponderomotive interaction with a leading quadratic 
dependence on position. 

Temporal harmonics in the lattice modulation drive 
atomic transitions at frequencies Hq = qUli, where q is 
an integer. The inset in Fig. illustrates how this is 
achieved. The EOM offset voltage is set to 14/2, where 
14 is the voltage difference between minimum and max¬ 
imum transmission. To modulate at the fundamental 
frequency Hi, one has V/ < 14/2 (dashed blue line). To 
access higher harmonics, V/ is increased (solid blue line). 
The resultant frequency upconversion is described by 


















2 



the incident intensity Lmc at the location of the atoms, 
scaled by Jdc (the intensity of the unmodulated high- 
power beam). Here, 77 is the power ratio between the 
modulated and unmodulated beams in the interferome¬ 
ter. A Fourier analysis of Eq. ([^ with 77 = 0.0077 (typical 
experimental value) leads to the Rabi frequency for the 
gth harmonic as a function of modulation strength 
(Fig. 1^. The Rabi frequency is scaled by I^c (in units of 
W/m^; -y/e, where e is the ratio of the return and inci¬ 
dent lattice-beam intensities (e = 0.09 experimentally); 
and D^f^. The unitless, normalized matrix elements 

for the spatial coupling, Z?" between Rydberg levels 
in the standing wave have been obtained in . For our 
lattice intensities and transitions, the Rabi frequencies 
are near 100 kHz. 

The Rabi frequencies exhibit a sinusoidal dependence 
on position in the optical lattice |4]; in Fig. we display 
the maximum Rabi frequency. The Rabi frequency for 
even-parity transitions (as in our case) is maximal for 
atoms at lattice intensity maxima and minima, where 
the leading spatial dependence of the ponderomotive po¬ 
tential is quadratic. Since the peak Rabi frequencies in 
Fig-i drop slowly as a function of q, unlike in typical 
nonlinear spectroscopy, we do not need to increase Idc 
to realize higher-order transitions, thereby avoiding in¬ 
creased light shifts. We utilize Fig. to determine the 
needed to achieve a high Rabi frequency at the har¬ 
monic order of interest. Atom-field interaction times 
are « 10 qLS. We target S ^ S transitions because 
they are insensitive to the MOT magnetic field (which 
is always on). Under traditional electric-dipole selection 
rules, these transitions would not be allowed in first-order 
perturbation theory. However, because we employ pon¬ 
deromotive spectroscopy, typical selection rules do not 
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FIG. 2. Expected scaled Rabi frequencies of the gth harmon¬ 
ics of the lattice modulation, as a function of V)i (units of 
W) i. See text for y-axis. 


apply [5]. 

In Fig. we demonstrate driving the ponderomo¬ 
tive transition 525'i/2 —t 53S'i72 via the third harmonic 
{q = 3) of the lattice modulation. For this transi¬ 
tion, Zlfff = 0.19, which is near the maximal value for 
nS {n + 1)5' transitions, = 0.215. The inset 

in Fig. [^a) shows that most spectral features can be 
reproduced by a semi-classical simulation |4]. The sim¬ 
ulation allows us to attribute the two outermost peaks 
to atoms that are anti-trapped (red-detuned peak) or 
trapped (blue-detuned peak) in the lattice. The simu¬ 
lation does not reproduce the sharp central peak in the 
experimental spectrum. This peak may be due to im¬ 
perfections in the optical lattice potential, which could 
arise from the lattice return beam encountering multi¬ 
ple interfaces that degrade the quality of the return focal 
spot. 

To determine the line center, we fit the smoothed, av¬ 
eraged spectral line in Fig. [^a) to a triple-Gaussian. To 
eliminate most of the trap-induced systematic shifts, we 
take the mean of the center locations of the outermost 
peaks. Assuming that the outermost peaks correspond 
to equal but opposite extremum light shifts (as would be 
the case in a perfect optical lattice), the mean provides a 
measurement of the line center with much-reduced light 
shift. See Table|T]for a comparison between our measured 
result, the calculated line center, and a two-photon mi¬ 
crowave spectroscopy reference measurement. 

In Fig. j^b) we show the dependence of the 535i/2 




FIG. 3. Third harmonic drive, (a) Population in 53S'i/2 as 
a function of lattice modulation frequency fli. Data are a 
smoothed average of 10 scans, 200 measurements each. Error 
bars, s.e.m. Red curve, triple-Gaussian ht. Vertical black line, 
line center. Black dashed line, center location of a two-photon 
microwave spectroscopy measurement (black curve). Inset, 
simulation results, (b) Peak 5351/2 population as a function 
of modulation strength V^. Each data point represents the 
peak height for an average of 6-10 scans with 200 measure¬ 
ments each. Vertical (horizontal) error bars, peak height (W) 
uncertainty. Red curve, proportional to the square of the 
q = 3 curve in Fig. 
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FIG. 4. Fifth harmonic drive, (a), (b) Population in 565'i/2 
as a function of lattice modulation frequency fli for Ryd¬ 
berg atoms prepared at lattice intensity (a) minima (b) max¬ 
ima. Data are a smoothed average of (a) 18 (b) 10 scans, 200 
measurements each. Error bars, s.e.m. Red curve, double- 
Gaussian fit. Vertical black line, line center. Insets, simula¬ 
tion results. 

excited-state population on modulation strength 
The behavior qualitatively agrees with the square of the 
Rabi frequency for the q = 3 curve in Fig. plot¬ 
ted for comparison. In particular, significant excited- 
state populations only occur beyond a threshold value of 
V^/14 « 0.6. Fig. j^b) reinforces that the spectrum in 
Fig. [^a) is due to the nonlinear q = 3 component of the 
lattice modulation. 

To reach higher-frequency transitions, it is desirable 
to drive transitions via even higher harmonics. In Fig.|^ 
we demonstrate the ponderomotive transition 54S'i/2 
56 «S'i/ 2, which has = 0.08, via the fifth harmonic 

{q = 5). The spectra in Figs. |^a), (b) are reproduced 
accurately by the simulations (insets). As before, the 
simulations indicate that the red-detuned (blue-detuned) 
peak represents atoms that are anti-trapped (trapped) 
in the lattice. In the figure, we also demonstrate the 
relation between Rydberg-atom position and light-shift 
polarity by either preparing the Rydberg atoms near an 
intensity minimum [Fig. Qa)], which results in mostly 
trapped atoms with positive light shifts; or maximum 
[Fig. ib)], which results in mostly anti-trapped atoms 
with negative light shifts [7]. 

To determine the line center, we fit the spectra in 
Figs. Qa),(b) to a double-Gaussian. We take the mean 
of the center locations of the outermost peaks, yield¬ 
ing a measurement in which the systematic light shifts 
mostly cancel. See Table [^ for results. Significantly, the 
q = 5 transition in Fig. j^has a transition frequency of 
about 94.4 GHz. Hence, we are approaching the sub- 
THz regime, which is important for improving ^vjv in 
precision frequency measurements. 

In spectroscopy, magic-wavelength lattices play an im¬ 
portant role because they allow for probing atoms in a 
light trap while avoiding systematic line-shifts due to 
the trap [9]. In a ponderomotive lattice, a magic condi- 
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FIG. 5. Magic condition, fundamental. Population in 70S'i/2 
as a function of lattice modulation frequency Di. Data are a 
smoothed average of 18 scans, 200 measurements each. Error 
bars, s.e.m. Red curve, single-Gaussian fit. Vertical black 
line, line center. Black dashed line, center location of a two- 
photon microwave spectroscopy measurement (black curve). 

tion occurs when a Rydberg atom becomes comparable 
in size to the lattice period, at which point the lattice- 
induced shifts of certain upper and lower Rydberg levels 
are nearly the same [10] and cancel. In a 1064-nm lat¬ 
tice, this phenomenon occurs for Rb nSi /2 atoms with 
lower- and upper-level principal quantum numbers sym¬ 
metric about 69.5. Furthermore, for principal quantum 
numbers between n = 66 — 73 the sign of the effective po¬ 
larizability of the atoms is reversed, indicating that Ryd¬ 
berg atoms are attracted to lattice-field maxima. This 
case differs from the more typical case where Rydberg 
atoms are attracted to field minima. 

In Fig. we show the results of driving the magic 
transition 69S'i/2 —t 70S'i/2, for which = 0.13, at 

the fundamental frequency {q = 1). The trap depths 
(and corresponding light shifts) for both states are about 
the same (2.2% of the free-electron ponderomotive trap 
depth). The magic condition results in an experimen¬ 
tal lineshape that is symmetric and has a narrow central 
feature. These characteristics are well-reproduced by the 
simulation (inset), which supports our interpretation of 
the signal. To determine the line center, we fit the cen¬ 
tral feature, expected to have nearly-zero systematic light 
shift, to a single Gaussian. As shown in Table |l| this 
measurement is determined to within 2 kHz (statistical 
uncertainty). 

For Fig. [^ we show a combination of both magic- 
condition and nonlinear ponderomotive spectroscopy, 
which has the greatest potential to improve jv. Here, 
we drive a magic ponderomotive transition 6851/2 
715 i/ 2, which has T^ggf = 0.08, via the third harmonic 
{q = 3). The trap depths of the 685i/2 and 715i/2 states 
are both approximately 2.0% of the free-electron pon¬ 
deromotive trap depth. In Fig. [^a), we prepare the 
Rydberg atoms at lattice intensity maxima and observe a 
narrow central peak on top of a broad, symmetric distri¬ 
bution; this central peak agrees very well with the result 
of a two-photon microwave spectroscopy measurement, 
also shown. In Fig. |^b), we prepare the atoms at lattice 
intensity minima and observe a narrow central dip on 
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FIG. 6 . Magic condition, third harmonic. Population in 
715 'i/ 2 as a function of lattice modulation frequency fli. Data 
are a smoothed average of (a) 12 (b) 5 scans, 200 measure¬ 
ments each. Error bars, s.e.m. Red curve, single-Gaussian fit. 
Vertical black line, line center, (a) Lattice not inverted. Black 
dashed line, center location of a two-photon microwave spec¬ 
troscopy measurement (black curve), (b) Lattice inverted. 

top of the broad distribution. In both cases in Fig.[^ the 
symmetry of the lineshapes agrees with our expectation 
for magic lattices. The spectral features are reproduced 
very well by the simulations (insets); the simulations also 
show that a central peak (dip) corresponds to trapped 
(anti-trapped) atoms. 

Because both Figs. [5][^ correspond to a magic condi¬ 
tion, we expect the respective frequency measurements to 
have nearly-zero light shift. We attribute the « 200 kHz 
discrepancy observed between Fig. |^a) and (b) to differ¬ 
ing residual electric fields, which may have caused differ¬ 
ent systematic DC shifts between the data sets (which 
were taken on different days). 

In Fig.[^(and to a lesser extent in Fig.|^, we note small 
oscillations near the central peaks. We attribute these 
to the effect of atoms running over the lattice potential 
wells, which are very shallow for the magic-condition lat¬ 
tices. For our atom temperature, atoms typically traverse 
multiple lattice potential wells during the modulation in¬ 
teraction time. Since the modulation-induced Rabi fre¬ 
quency flips sign twice per lattice period, a rotary-echo¬ 
like effect m occurs. This causes periodic structures 
when varying the detuning. 

An important finding in both Figs. [5][^ is that the 
central peaks in the lattice-modulation spectra and the 
peaks in the two-photon reference spectra are at Fourier- 
limited resolution. Therefore, precision measurements 
made via magic-condition ponderomotive spectroscopy 
and traditional microwave spectroscopy can have simi¬ 
lar spectral resolution. However, ponderomotive spec¬ 
troscopy allows us to access a wide variety of typically- 
forbidden transitions at high frequencies. 


_gHi (GHz) i^caic (GHz) 2 t/ 2 p (GHz) 

52Si/2 5351/2 54.693423(15) 54.693577(5) 54.693556(6) 

q = 3 

54 Si /2 ^ 56Si/ 2 a. 94.370915(15) 94.371060(8) 
q = 5 b. 94.370955(15) 94.371060(8) 

69Si/2 70Si/2 22.509226(2) 22.509227(1) 22.509226(2) 

q = 1, magic 

68 S 1/2 7 IS 1/2 a. 67.588473(2) 67.589048(4) 67.588478(1) 

g = 3, magic b. 67.589046(9) 67.589048(4) 

TABLE 1. Summary of results. All measurement uncertain¬ 
ties, statistical. See text for details. 

Frequency-measurement results are summarized in Ta- 
ble|i| where we list transition frequencies (gHi) measured 
in the modulated lattice, expected transition frequencies 
(t'caic) calculated using quantum defects from [H], and 
reference measurements ( 2 ^' 2 p) obtained via two-photon 
microwave spectroscopy in the absence of a lattice. Over¬ 
all, agreement is quite satisfactory. Statistically signifi¬ 
cant shifts in the measured results from the calculated 
values are negative, indicating possible Stark shifts due 
to residual DC electric fields. In the magic-condition 
cases, the statistical uncertainty in gHi is of the same 
order of magnitude as the uncertainty in 2r'2p. With im¬ 
provement of stray electric field control, the modulated 
lattice should be suitable to measure atomic transition 
frequencies that are difhcult or forbidden in microwave 
spectroscopy. 

In conclusion, we have performed ponderomotive spec¬ 
troscopy of Rydberg atoms by employing higher harmon¬ 
ics of the lattice modulation as well as magic conditions 
of the lattice. We have reached transitions near the sub- 
THz regime, with the potential to exceed that limit, using 
lower-frequency microwave sources. Using magic lattices, 
we have demonstrated Fourier-limited spectral lines with 
a small Az/, as well as reduced light shifts. Being able 
to access sub-THz atomic transitions with low Av and 
free from light shifts will improve the fractional resolu¬ 
tion (^Avjv) of transition frequency measurements. Ap¬ 
plications include precision measurement of atomic char¬ 
acteristics m and physical constants (e.g. the Rydberg 
constant [HI, which could lead to verification of the pro¬ 
ton size m)- 
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